Tracking fast small color dipoles through strong gluon fields at the LHC 
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We argue that the process 7 + A — > J/tp + " gap" + X at large momentum transfer q 2 provides a 
quick and effective way to test onset of a novel perturbative QCD regime of strong absorption for 
the interaction of small dipoles at the collider energies. We find that already the first heavy ion run 
at the LHC will allow to study this reaction with sufficient statistics via ultraperipheral collisions 
hence probing the interaction of qq dipoles of sizes ~ 0.2fm with nuclear media down to x ~ 10~ 5 . 
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Soon after J/ip was discovered, the J/ip photoproduc- 
tion experiments on nuclear targets have established that 
nuclei are practically transparent to the J/ip's produced 
at photon energies in the range ~ 20 GeV -j- 120GeV 
[TJ The absorptive cross section, cr^f N , was found 
to be close to ~ 4mb that is much smaller than the 
cross section of interaction of ordinary mesons ~ 25 mb. 
The observed transparency is natural within the Low- 
Nussinov model of two gluon exchange where the cross 
section of hadron interaction with a small color singlet 
dipole quark-antiquark configuration in the photon wave 
function is proportional to the square of the transverse 
size of the color dipole 011]. Note, that the average size 
of cc configurations involved in photoproduction of J/ip 
is significantly smaller than the J/ip size. Such suppres- 
sion of interactions of small dipoles is well known effect in 
electrodynamics - for example, a muonium can propagate 
through the media much easier than a positronium. 

Within the leading InQ 2 , ln(l/a;) approximations of 
perturbative QCD one expects (in difference from the 
Low - Nussinov model) that the cross section of the in- 
teraction of small dipoles with hadrons should increase 
rapidly with increase of invariant dipole-hadron energy 
WryN — y/s due to the growth of the gluon fields in hadron 
targets at small x oc s _1 : 

a dlp . T (x, d) = — F 2 d 2 a s (\/d 2 )xG T (x, \/d 2 ), (1) 

where F 2 — 3(4/3) is the Casimir operator for the two- 
gluon (qq) dipole, and A ~ 4 — 9. For a dipole of a size 
~ 0.25 fm relevant for production of J/ip Eq.Q corre- 
sponds to energy dependence cx s ' 2 and describes well 
behaviour of the both exclusive electro production of vec- 
tor mesons and the inclusive cross section of deep in- 
elastic electron-proton scattering scattering observed at 
HERA, for a review and references see [5]. 

A naive extrapolation of the observed pattern to LHC 
energies indicates that the strength of this interaction 
may reach values comparable to that experienced by light 
hadrons, leading to a new regime of strong interaction 
physics at the LHC characterized by a strong absorption 
of small color dipoles by the media. On the other hand 



it is evident that to avoid conflict with probability con- 
servation starting from some energies such rapid increase 
of cross section should be tamed. 

So, the question is whether it will be possible at LHC 
to observe this new perturbative QCD regime when the 
coupling constant is small but the interaction is strong. 
In practice, it is very difficult to devise a high energy 
probe for virtualities of few GeV 2 for the hadron col- 
liders especially for the high energy strong interactions 
involving nuclei where gluon densities per unit area are 
higher and where new high gluon density physics should 
be enhanced. Such a problem is absent for electron - ion 
colliders but these colliders are far in the future. 

An alternative which we discuss here is to use ultrape- 
ripheral collisions (UPC) of ions at the LHC in which 
one of the nuclei serves as a source of quasireal pho- 
tons and another one as a target. The recently pub- 
lished study [B] demonstrates that it is feasible to se- 
lect UPC at the LHC and that the rates for many pro- 
cesses of the dipole - nucleus interactions are high enough. 
This includes the process of coherent photoproduction of 
-//^d- However, this process could only be effectively 
studied up to relatively small energies W 7 jv ~ 130 GeV 
due to the inability to separate contributions due to the 
lower and higher energy photons emitted by two colliding 
ions. Here we suggest a strategy which avoids the above 
mentioned shortcoming of the coherent J/ip production. 
It is based on the study of the large momentum transfer 
-t= q 2 = (P'f-pj/4,) 2 process: j + A — > J/ip+ gap + X. 
In addition to the theoretical advantages which we will 
explain below it also has some appealing observational 
features. Observation of J /tp and hadrons allows to de- 
termine unambiguously which of the nuclei emitted the 
photon. As a result it is possible to observe the process 
up to W 1 -n ~ ITeV. Besides, acceptance of the all three 
LHC detectors which plan to study heavy ion collisions 
is sufficiently large for the discussed kinematics. 

At first, let us briefly consider the main features of the 
elementary reaction 7 + p — > J ftp + gap + X in a case of 
a proton target which was extensively studied at HERA. 
This process belongs to a class of reactions introduced 
in [5] [9] where selection of large q 2 ^> ensures that 
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FIG. 1: The typical perturbative QCD diagram for 
J/ip production with gap. 



transverse momenta remain large in all rungs of the gluon 
ladder and therefore two gluons are attached to one par- 
ton in the target, see Fig. [T] As a result the cross section 
can be written in a factorized form as a product of the 
gluon density in the target (the kinematics can be chosen 
so as the correction due to scattering off the quarks which 
we will not write below explicitly will be small) and the 
elastic cross section of scattering of a small dipole off a 
gluon: 



-yT^J/^X 

dq 2 dx 



da 



dq 2 



9T(x,q 2 ), 



(2) 



is the gluon distribution in the target, 



where gr(x, q 2 ) 

x = q 2 /(M| — m 2 N + q 2 ) and Mx is the invariant mass of 
the produced hadron system. Though it is not practical 
to measure M x directly, it can be expressed with a good 
accuracy through Ay, the rapidity interval occupied by 
the system X as 



Ay = l^y/q 2 /xm N ). 



(3) 



The HERA data for the large q 2 and rapidity gap pho- 
toproduction of J/ip off the proton target are consistent 
with the expectations of the models based on the factor- 
ization approximation in Eq.Q and QCD inspired mod- 
els for the 7(7 — > J/ipg amplitude, see Refs. [ID], [TT] for 
recent data analyses and references. 

Now let us analyse the effects arising in a case of the 
nuclear target in which the color singlet dipole produced 
in the photon interaction with one of the nuclear nucleons 
passes through nuclear medium before to form the J / tp . 
Firstly, we need to estimate the average dipole size, d, 
which enters in the elementary 75 — > gJ/ip amplitude. 
The overlapping integral between photon and charmo- 
nium wave functions is determined by c-quark momenta 
k t ~ m c . For a two body cc system with Coulomb inter- 
ior q 2 <C m 2 , which is 



action one finds d 
similar to the estimate of d » .25 fm [12] for the process 
of forward J/ip coherent photoproduction with a charmo- 
nium realistic wave function. For large q 2 in the Coulomb 
atom approximation we find d 2 (q 2 )/d 2 (0) w Am 2 /q 2 for 
the dominant amplitude when gluons are attached to the 
same c-quark. Hence we will use in our analysis the in- 
terpolation formula 



with do = .25fm,m c = 1.5 GeV. 

We choose the kinematics where x > x s h = 0.01, 
hence, the shadowing and antishadowing effects in the 
nucleus parton density in the process of hard interaction 
are small, gA(x,q 2 ) « Ag p (x,q 2 ), and Eq.Q leads to a 
cross section which is approximately linear in A. 

Now we demonstrate that diagrams where the ladder 
is attached to gluons belonging to different nucleons of 
the nucleus give a small contribution for realistic nuclei in 
spite of been enhanced by a combinatorial factor cx A 1 / 3, . 
For simplicity we consider scattering at central impact 
parameters for the rapidity interval given byEq.jjj). In 
this case one of the partons should carry (3i light cone 
fraction of the nucleon momentum, which is comparable 
to x, while the second one is allowed to have any /?2 > fti- 
To obtain an upper limit for the ratio, R2, of contribu- 
tions of these diagrams to that described by diagram in 
FigjTJwe assume that momentum transfers in the interac- 
tions with two partons of the nucleus are q\ ~ q\ ~ g 2 /4. 
Emission of two partons with p t = \fq 2 /2 to the same 
angular interval as for the single parton with p t = yfq 2 
requires (3\ = x/2 and P2 > Pi- Hence in this approxi- 
mation i?2 is expressed through the number of gluons in 
the cylinder of radius d/2 with f3\ — x/2, 02 > 01 as (cf. 
Eq.(l)) : 
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where T(b) = J_ pA(b,z)dz is the standard thick- 
ness function normalized to A. A factor of two in the 
nominator is due to the presence of two different at- 
tachments of gluons to partons of the target. For the 
J/ip production, using Eq.Q we find, for example, for 
q 2 = 4GeV 2 ,x = x sh , R 2 < 5%. The factor R 2 stays 
small at larger i due to decrease of parton densities with 
x. The same conclusion is valid for x > 0.1 where contri- 
bution of scattering off quarks becomes noticeable. Note 
also that the momentum balance in the case of interac- 
tion with two different partons is strongly different than 
in the case of interaction with a single parton leading to 
possibility to suppress this contribution using kinematic 
cuts. 

Hence we conclude that, in the discussed kinematics, 
significant deviations from the linear A-dependence can 
arise only due to possible inelastic interactions of the 
dipole with other nucleons of the nucleus. The change of 
the probability of the gap survival with A: 
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d 2 (q 2 )/d 2 (0) w (l + (z 2 /4m 2 )- 



(4) 



which is formally a higher twist effect in d 2 has a physi- 
cal meaning of the probability for a small dipole to pass 
through the media without inelastic interactions for the 
energy W lN . 
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This effect can be expressed through the profile func- 
tion for the dipole - nucleus scattering, T d i P .A{x, d, b) , 
which is the Fourier conjugate of the elastic dipole - A 
amplitude. It is normalized so that cr tot (dip — A)(x, d) = 
2 J d 2 bT 'dip,A(x,d,b). The range of gluon x probed in 



this case is of the order x = ^jm/H^jv ( an d some- 



what smaller if one uses the charmonium model for the 
J/ipwave function like in |12j). In the dynamics driven 



by inelastic interactions 



< 1. Application of 



T(x, d, b) 

S-channcl unitarity (essentially the probability conserva- 
tion, cf. [13] ) allows one to demonstrate that the prob- 
ability for the dipole not to interact inelastically is equal 



to 



i - T(x, d, b) 



leading to 



A eff = / d 2 bT(b) 1-T(x,d,b) 



(7) 



Here we neglect fluctuations in the size of the dipole 
which is a good approximation for the regime of moder- 
ate absorption where average interaction strength enters 
into the answer. In the case of large absorption the fil- 
tering effect takes place leading to enhancement of the 
contribution of small dipoles. A more detailed treatment 
will be given elsewhere. 

Choice of kinematics with x > x s h results in domi- 
nance of hard interaction at small impact parameters. 
Thus, using a heavy nucleus as a target one can probe 
propagation of a small dipole through ~ 10 fm of nuclear 

matter and determine 1 — T(x, d, b) . 

To estimate the suppression effect as given by EqjT] 
we use two popular complementary models for the in- 
teraction of a small size dipole with the matter. One is 
the eikonal model where the small size dipole interacts 
via multiple rescatterings off nucleons with the strength 
given by the dipole - nucleon total cross section. Devi- 
ations of the dipole - nucleus interaction from oc A is 
a higher twist effect since the interactions with n > 2 
- nucleons is oc d 2n . Second model is the leading twist 
shadowing model which includes only two gluon attach- 
ments to the dipole. In this case deviations from the 
linear regime in A are due to soft interactions of the two 
gluons with the nucleus. 

In the eikonal model, neglecting fluctuations of the cc 
transverse size we obtain: 

r(x, d, b) = 1 - exp(-a dip - N (x, d)T(b)/2), (8) 



where <Jdi P -N(x,d) is given by Eq.pJ). Since for heavy 
nuclei T(0) re 2fm" 2 , Pf v w exp[^a dip ^ N (x, d)T(O)) 
becomes small already for a ~ 5 mb which corresponds 
to x ~ 1CP 3 . Hence in this model a large suppression 
effect is expected which grows with W^m and, for fixed 
W-yN, decreases with increase of q 2 , see Fig. 2 (the curves 



for q 



50GeV aim to illustrate the trend of the t- 



dependence of -Pf ap , the actual measurement for this 
range of t will require a long running time) . 



An alternative model is the leading twist approxima- 
tion over parameter Aq CD / (Am 2 +q 2 ) for the dipole scat- 
tering off the nucleus which was used for the description 
of coherent J/ip production. Contrary to the eikonal ap- 
proximation this approach accounts for essential nuclear 
modification of the nuclear parton distributions at small 
x. Since in the leading twist Eq.([lJ describes inelastic 
dipole - nucleus cross section, the probability for a dipole 
of the size d to pass through the nucleus without inelastic 
interactions is 



PT P =J I d 2 bT(b)[l 



-N 



(x,d) 



g A (x,Q 2 ,b) ^ 
gN(x,Q 2 ) 



(9) 



where Q 2 = X/d 2 and gA(x, Q 2 , b) is the gluon 
density of the nucleus in impact parameter space 
(/ gA(x,Q 2 ,b)d 2 b = gA(x,Q 2 )). In the kinematic range 
x > 3 • 1CP 3 where shadowing effects are still small one 
can unambiguously calculate the shadowing correction as 
a function of b through the diffractive gluon parton dis- 
tribution function (pdf), gdiff{x,xjp,Q 2 ) which is mea- 
sured in hard processes at HERA. Higher order rescat- 
terings could be estimated by introducing a e ff(x, Q 2 ) — 

J x dxpg d rff{x 1 X]p,Q 2 )/g N (x,Q 2 ), for details see [H]. 

The very small x and low virtuality diffractive pdf's 
one has to use for such an analysis are not reliable as 
they involve extrapolations from larger Q 2 and x. A 
straightforward application of the data leads to a very 
strong shadowing of gA and hence to a small absorption, 
sec Fig{2] (dotted line). However it is very difficult to 
envision leading twist (LT) dynamics where partons of 
nucleons at a given impact parameter, 6, would screen 
the nuclear pdf below the maximal value of generalized 
gluon density gN(x, Q 2 , p) of one nucleon at this b (in the 
Glauber model for the nucleon-deuteron interaction this 
condition corresponds to a(hD) > a(hN)). For the limit 
of large A this implies a condition 



g A (x,Q 2 ,b = 0) >g N (x,Q 2 ,p=0). 



(10) 



An effective way to implement this condition is to use 
the eikonal expression for screening of Qa{x, Q 2 , b): 



g A (x,Q 2 ,b) = 2 
gN(x,Q 2 ) <T eff (x,Q 2 ) 



[1 -exp(- 



-T(b)a eff (x,Q 2 ), 



which leads to gA(x,Q 2 ,b)/gN(x,Q 2 ) < 2/a e ff(x,Q 2 ) 
and find maximal value of a e f f(x,Q 2 ) which satisfies 
Eq 10 for large T(b). The t-dependence of the gluon 
generalized parton distribution g^(x,t) is measured at 
HERA in the exclusive vector meson production. For 
x - 10~ 4 exponential fits (exp(Bi)) find B - 5GeV" 2 . 
Using this fit we find a^ff = AnB - 27 mb. This ge- 
ometrically constrained LT (gcLT) model leads to the 
solid curve in Fig(2] One can see that the magnitude of 
the suppression predicted by such gcLT model is rather 
close to that obtained in the eikonal model, the effect of 
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FIG. 2: The rapidity gap survival probability as a function of 
W jN for q 2 = 4 GeV 2 and q 2 = 50 GeV 2 

To determine for what kinematic range at LHC the 
measurements of cross section of the ultraperipheral pro- 
cess A + A — > A+ J/ip + gap + X are feasible we used 
the standard expression for the photon flux generated 
by the relativistic ions. To calculate the cross section of 
7P — > J/tp + gap + X we adopted the QCD motivated 
parametrization of [10 which fits all relevant HERA data 
using CTEQ6L parton distributions [15]. P^ op was cal- 
culated within the models described above. We imposed 
the cut: x m in > 0.01 corresponding to Mx < 10y/qf*. 

Since it is possible to determine experimentally which 
colliding ion serves as the source of photons we don't 
account here for the symmetrical contribution from an- 
other ion which increases the counting rate by a factor 
of two. The results presented in Fig (3] indicate that for 
the planned integrated luminosity for one month of run- 
ning per year 4.2 x 10 5 mb _1 , it will be possible to mea- 
sure P Q j^ v in a wide range of q 2 up to W^n ~ 1 TeV. 
A much larger range of q 2 will be feasible for the pro- 
cess 7 + A — ► p + gap + X both due to a larger el- 
ementary cross section and absence of the suppression 
due to a small branching ratio for J/ip — ► . Mea- 
surements with p's will also allow to determine up to 
what q 2 the end point contribution which corresponds to 
very small P^ ap cx A~ 2 ^ 3 gives a noticeable contribution 
and at what q 2 the leading twist contribution leads to 

In conclusion, we find that a study of the J/ipm UPC 
with rapidity gaps will allow to determine the pattern 
of absorption for the interaction of small dipoles with 
0.15/m < d < 0.25 fm up to W lN ~ ITeV which corre- 
sponds to interaction with gluon fields with x down to 
10~ 6 for virtualities of ~ 4 4- 10 GeV 2 . The ability to 
vary both q 2 and W^n will allow to separate the regions 
where opacity should be large and where the high energy 
color transparency will hold. Extension of these measure- 
ments to smaller q 2 will help to get better understanding 



of soft diffractive dissociation induced by a small dipole. 
The discussed J /ip measurements will also serve as an im- 
portant benchmark for the studies of J /if) production in 
the heavy ion collisions at the LHC as well help to un- 
derstand the possible role of small dipoles in generating 
strongly fluctuating pedestal in the inelastic pp collisions. 
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